Antiferromagnetic metal to heavy-fermion metal quantum phase transition in the 
Kondo lattice model: A strong coupling approach 
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We study the quantum phase transition from an antiferromagnetic metal to a heavy fermion metal 
\^ , in the Kondo lattice model. Based on the strong coupling approach we first diagonalize the Kondo 

■ coupling term. Since this strong coupling approach makes the resulting Kondo term relevant, the 
Kondo hybridization persists even in the antiferromagnetic metal, indicating that fluctuations of 

i Kondo singlets are not critical in the phase transition. We find that the quantum transition in our 

j i ■ strong coupling approach results from softening of antiferromagnetic spin fluctuations of localized 

' spins, driven by the Kondo interaction. Thus, the volume change of Fermi surface becomes contin- 

uous across the transition. Using the boson representation of the localized spin fti = \z\ a T tT<T i z ia i 
with the spin-fractionalized excitation z, CT , we derive an effective U(l) gauge Lagrangian in terms 
of renormalized conduction electrons and fractionalized local-spin excitations interacting via U(l) 
gauge fluctuations, where the renormalized conduction electrons are given by composites of the con- 
duction electrons and fractionalized spin excitations. Performing a mean field analysis based on 
i—i ■ this effective Kondo action, we find a mean field phase diagram as a function of Jk/D with var- 

ious densities of conduction electrons, where Jk is the Kondo coupling strength and D the half 
j_( , bandwidth of conduction electrons. The phase diagram shows a quantum transition, resulting from 

■ condensation of the spin-fractionalized bosons, from an antiferromagnetic metal to a heavy fermion 
metal away from half filling. We show that beyond the mean field approximation our critical field 
theory characterized by the dynamic critical exponent z — 2 can explain the observed non-Fermi 

S liquid physics such as the specific heat coefficient 7 = C v /T ~ — InT near the quantum critical 
point. Furthermore, we argue that if our scenario is applicable, there can exist a narrow region of 
>-^j ' an anomalous metallic phase with spin gap near the quantum critical point. 

c : 

PACS numbers: 71.10.-w, 71.10.Hf, 71.27.+a, 75.30.Mb 
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^ ■ I. INTRODUCTION 

> 

Nature of non-Fermi liquid physics near quantum critical points in heavy fermion compounds is one of the central 
interests in modern condensed matter physics. [1] A standard theoretical framework is the Hertz-Moriya-Millis (HMM) 
theory in the Landau- Ginzburg- Wilson framework. [2] Using a Hubbard-Stratonovich (HS) transformation for an 
appropriate interactional channel, a local order parameter can be introduced. Integrating out fermion degrees of 
freedom and expanding the resulting logarithmic action based on the noninteracting itinerant fermion ensemble, one 
' obtains an effective action of the local order parameter with dissipation that results from gapless electrons near the 
"jS ■ Fermi surface. Based on this order parameter action, a self-consistent mean field analysis and a renormalization group 
study can be performed to find non- Fermi liquid physics near the quantum critical point, characterized by critical 
fluctuations of the local order parameter. [3] 

However, since the expansion in deriving HMM theory from the Fermi liquid action is basically a weak coupling 
approach, the effective theory will break down in a strong coupling limit. In addition, the presence of gapless 
electrons can cause nonlocal interactions between order parameters, making a conventional treatment unreliable in a 
local effective action. [4] Furthermore, HMM theory is not fully self-consistent because feedback effects to the fermion 
• • ■ degrees of freedom by order parameter fluctuations are not taken into account. Actually, it is now believed that the 
critical field theory of order parameter fluctuations, that is, HMM theory cannot explain the observed non-Fermi 
liquid physics in thermal and electrical properties in the heavy fermion metals. [1] 

In this study, using a strong coupling approach, we derive an effective action in terms of both bosonic and fermionic 
excitations associated with order parameter fluctuations and gapless electrons, respectively. The weak coupling 
approach solves the kinetic energy term first, and treats the interaction term perturbatively, implying that the theory 
is based on the noninteracting itinerant fermion ensemble. On the other hand, in the atomic limit of electrons the 
strong coupling approach solves the interaction term first and then, treats the kinetic energy term perturbatively. In 
the present strong coupling approach we take into account critical order parameter fluctuations and gapless electron 
excitations near the Fermi surface on an equal footing. In other words, we incorporate both bosonic and fermionic 
excitations in the effective action without integrating out the gapless fermion degrees of freedom. Hence, our strong 
coupling approach goes beyond the conventional treatment for the quantum phase transition in interacting itinerant 
electrons. 

We start from an antiferromagnetic phase of localized spins in the Kondo lattice model, where the Kondo interaction 
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term describes the coupling between the localized spins Hi and conduction spins s, = \(y ia T aa iCi (J i . In the strong 
coupling approach the Kondo interaction term can be diagonalized by using the CP 1 representation Hi • r = UiTzU} 

( -z f \ 

with an SU(2) matrix Ui — I f ) ^ or tne l° cauze d spins and introducing renormalized electrons Via as 

composites of the conduction electrons and spin-fractionalized bosons, tp icr = U^ aa , Ci a > . [5-7] Then, the Kondo lattice 
model can be rewritten in terms of the spin-fractionalized bosons Zi a and renormalized conduction fermions tpi a . 
Integrating out the fermion fields ipia and performing the gradient expansion in the resulting logarithmic action, 
this conventional strong coupling theory leads to an effective action of the Zi„ bosons, which is well known in the 
context of the nonlinear a model. [5] In the present study, however, we take into account both the bosons and fermions 
simultaneously. 

The boson representation for the localized spin instead of its fermion representation is necessary for explaining 
the quantum phase transition involving an antiferromagnetic phase. The U(l) slave-boson representation has been 
used conventionally for the study of the heavy fermion phase in the Kondo lattice problem. If the localized spins are 
represented by fermions and an order parameter (slave-boson) corresponding to the Kondo hybridization is introduced, 
a critical theory can be obtained in terms of both order parameter fluctuations and fermion excitations by integrating 
out the conduction electrons. [8, 9] Fermions and bosons then interact via long range interactions mediated by slave- 
boson U(l) gauge fields. This critical theory successfully explains the non-Fermi liquid physics such as the specific 
heat coefficient 7 = C v /T <~ — InT near the quantum critical point. However, it is difficult to explain the observed 
antiferromagnetic long range order that begins to appear at the transition point where the heavy fermion phase 
disappears. Furthermore, the order parameter in the slave-boson representation is "hidden"; no symmetry breaking 
associated with lattice translations or spin rotations is involved. Thus, the transition driven by condensation of the 
hidden order parameter may not be physical. 



II. EFFECTIVE ACTION FOR THE KONDO LATTICE MODEL 



In the Hamiltonian of the Kondo lattice model 

HkLM = ~<yi Ao C 3° + JK^Si- C^Taa'Cia 1 , (1) 
ija iaa' 

the first term describes dynamics of conduction electrons Ci a and the second term represents antiferromagnetic ex- 
change couplings between conduction electrons and localized spins 5^ , where t is the hopping integral of the conduction 
electrons and Jk is the Kondo coupling strength. Using the coherent state representation for the conduction electrons 
and localized spins, [10] the partition function of the Kondo lattice model in the path-integral representation can be 
given by 

Z = J Dc ia Dn i c^p\iS'^ J dr J durii ■ x 

- I dT (Yl c ^ dr ~ ^ Ci ° ~ 1 c i° c 3° + JkS Y Hl ' i^"'^') ' ( 2 ) 



i] a 



where we use Si = Srii with S =1/2 and |nj| = 1, and \i is the chemical potential. The first term corresponds to Berry 
phase that comes from the path- integral quantization of spin with an additional parameter u in a unit sphere. [10] 
In order to study the antiferromagnetic phase with collinear ordering of localized spins, we set 

n t -» {-l) l n t . (3) 

Then, Eq. (2) reads 

Z = f Dci a Dniexp[iSj2(-iy [ ^ j\un t ■ x ^) 

-/ dT (j2 C l^ - V^a -tJ2 C ia C j* + ■fr<S$^(-l) i ct r (fr i • rWc iCT <)] . (4) 
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Based on this partition function we investigate two closely-related interesting problems; how the antiferromagnetic 
ordering of local spins is affected by dynamics of conduction electrons as varying the Kondo coupling , and how the 
dynamics of the conduction electrons is influenced by the change of antiferromagnetic fluctuations of the local spins. 

We apply the strong coupling approach to Eq. (4) where we first solve the Kondo coupling term. Although the 
HMM theory results in a firm and successive quantum theory for quantum phase transitions, [2] the naive perturbation 
for the Kondo coupling term based on the itinerant fermion ensemble do not work well in the case of strong couplings. 
Using the identity 



■f=U t T 3 u}, Ui = I Zn ! a > (5) 




where Zi a is a boson field carrying a spin a in the SU(2) matrix field Ui and satisfies the unimodular constraint 
J2a \ z iv\ 2 = 1' an d performing the gauge transformation 

i) ia = ut aa ,c ial , (6) 

the Kondo coupling term, J K S J2iaa> ( _ l) lc L(™j ' T) aa >c iG > , can be represented as JKSJ2iaa>(~ 1 Y' l l>h T vcr>iPi>?'- As a 
result, the two-body Kondo interaction term is represented by a one body term. We call Zi„ and ipi a bosonic spinon 
and fermionic chargon, respectively. 

In the strong coupling approach an antiferromagnetic spin fluctuation Hi carrying spin quantum number 1 fraction- 
alizes into bosonic spinons Zi„ with spin 1/2, which seems to occur through the screening of conduction electrons due 

fa J V -^i c »T + z n c n 

where t/^t field represents the usual Kondo hybridization, and tpi^ field the polarization of the bosonic spinon and 
the conduction electron. The fermions ipia- can be considered to express Kondo resonances. Another way to describe 
this fractionalization is that the conduction electrons fractionalize into the bosonic spinons Ui aa ' and the fermionic 
chargons ^W, i.e., c ia = Ui aa '^i a i . The resulting partition function is obtained in terms of new field variables -0i and 



to the strong Kondo interaction. The components of ipi a field are given by ipi„ = 



Z = J Dip la DU taal 5(U} IJiy „U t(Tll(jl - <W)exp 



ij icrcr' 

where Sb is the Berry phase action. Note that the integration measure J DciuDfii in Eq. (4) is changed into 
/ Dip i(T DUi aa 'd(Uj rTrJ „Ui a " (7 ' — 5 arr >) in Eq. (7). Note that the chargons (renormalizcd conduction electrons) and 
spinons (fractionalized local spins) are now coupled in the kinetic energy term instead of the Kondo coupling between 
the conduction electrons and localized spins in Eq. (4). 

A standard way to treat this nontrivial kinetic energy term is to integrate out the chargon fields, 



DUi aal 5(U ] iaa „U i<T » a > - cW)exp 



-S 



B 



+trln 



{[d T - ,i]<w + J K s{-i) l Tl a , + [u}d T Ui] aa , - ^t/Ltw}] , (8) 



where tij = t is the nearest neighbor hopping. An effective action of the spinons can be obtained by expanding the 
logarithmic term for the bosonic spinons. One important difference from the HMM theory is that the expansion 
parameter is t/Jjc instead of Jx/t. However this formulation has a serious defect. Metallic physics of the conduction 
electrons is not introduced since this treatment is valid in the atomic limit. Actually expanding the logarithmic term 
in the expansion parameter tj 3k, the resulting effective action is known to be the 0(3) nonlinear a model appropriate 
to an insulating antiferromagnet.[3] Because an additional Berry phase term appears in the effective a model, the two 
Berry phase terms cancel each other and the contribution of Berry phase vanishes. [6] In the following, although we 
develop a formulation different from the above standard approach, we can also exclude the Berry phase term. 

An important issue of this study is how to introduce physics of the Fermi surface of the conduction electrons in the 
strong coupling approach. One possible route is to decouple the "interacting" kinetic energy into the conventional 
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" nonintcracting" one via the HS transformation. Unfortunately, there is a difficulty in performing the HS transfor- 
mation for the time-derivative term in Eq. (7). We consider discrete-time steps and rewrite the partition function 

as 

itt' ijr 

+^E(-i)V£^ -jE^c}], (9) 

ir ir 

where J is an energy scale associated with a time step J = 1/At. This discrete-time expression can be reduced to 
the original one of Eq. (7) in the limit of At — > 0.[11] 

Performing the HS transformation for the "hopping" terms in Eq. (9), we obtain the following expression 

irr' ijr 

= f DF™'DE™' exp[- Y,{E\^Ffrr> + H.c U^U^' F™, E^,^t H - c ) 

J irr' 
ijr 

= f DF^DE^exp[-J2tr{El T ,F iTT , + H.c. - XJ iT ,F iTT ,u\ r - ^\ r E\ TT ^ lT , - H.c.} 

irr' 

~ 'J E tr { E hrFijr + H.C. - U jT F ijT Uj T - ^Mjr^jr ~ H.C.}] , (10) 



where E irT > , F iTT > and E^ T , F^ T are HS matrix fields associated with hopping of tp iT fermions and z iT bosons in time 
and space, respectively. 

We make an ansatz for the hopping matrix fields 

E iTT , = E t e ia ^' T \ E ijT ee E r e ia ^ r \ 

F lTT , ee Fte**"' 7 *, F ijT = F r e ia ^ T3 , (11) 

where Et, Ft and E r , F r are longitudinal modes (amplitudes) of the hopping parameters, and ai TT ' , aijr their transverse 
modes (phase fluctuations) which are time and spatial components of U(l) gauge fields. In fact, the transverse modes 
should be represented by SU(2) gauge fields generally but we limit our discussion in the U(l) case for simplicity. 
The hopping parameters will be determined self-consistently. Inserting Eq. (11) into Eq. (10), and using the explicit 
representation [Eq. (5)] of the SU(2) matrix [/,, we obtain an effective U(l) gauge theory of the Kondo lattice model 

Z = j Dz lTa Di, lTa Da^D\ lT DE t DF t DE r DF r e^[-{^^E t F t + - j Y,Y. E - rF - r 

tt' i t ij 

E E ipLe iaai -'Ar'« -jE r e E 4r^° aii ^ir„ + ^ E B-i^L^w 

tt' ia t ija r ia 

- -j E E tL**™ - Ft E E *»* -j^EE4 

r ia tt' ia r ija 

+*EE^(Ei^i 2 -i)}]. ( 12 ) 

t i o- 

where the T3 matrix is replaced with a = ±. is a Lagrange multiplier field to impose the unimodular constraint. 

The last step in deriving an effective U(l) gauge theory of the Kondo lattice model is to perform the limit of 
At — > 0. Ignoring the time component ai TT ' of the U(l) gauge field for the time being, the time-derivative term of 
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the bosonic spinons becomes in the tight-binding approximation for the discrete time 

a, 
1 



--2F^^(l-i(^) 2 )4(0„), iCT (n„) 

= [ P dTj2(^\drZ itT \ 2 -2JF t \z ta \ 2 ). (13) 
•to lff w 7 

The last term can be absorbed into the Lagrange multiplier term, thus has no physical effects. The derivation of 
Eq. (13) is based on the relativistic invariance for the bosonic spinons. This is reasonable because spin excitations 
in the antifcrromagntic phase have the uj n ~ k dispersion, exhibiting the relativistic invariance. On the other hand, 
the relativistic assumption for the bosonic spinons is not appropriate for the fermionic chargons because the chargons 
have Fermi surface. In this case it is natural to perform the At — > limit naively. Then the time-derivative term for 
the fermionic chargons is given by 

-£*EE tir^<° = -£*EE *L 0- - ^ff at) 

tt' ia t ia 

= f dT E (^lAVw - JEk^J) , (14) 

where Ar = 1/J is used. The last term is also absorbed into the chemical potential term. 

Based on the above discussion, we find the effective U(l) gauge action for the Kondo lattice model 

Z = j DziaDipurDa^DXiDEtDFtDErDFre- 8 '", 

&eff = So + S,p + S z , 

So = f d T (jJ2EtF t +tJ2 E r F r ) , 

i ij 

S,= dT[E t J2 4*{dr - ioa lT )^ la - tE r i>le iaaii i>j* + 4S^(-l)V4fc 

^ ia ija ia 

-me^^)' 

ia 

S z - J dr(^^|(9 r - ia lT )z la \ 2 ~tF r J24*^z ja +iJ2^(J2 ~ ^) ■ ( 15 ) 

^ icr ija i a 

Here, both the conduction electrons and localized spins are taken into account on an equal footing in the strong 
coupling regime, and the bosonic effective action S z associated with order parameter fluctuations is derived without 
integrating out the gapless conduction electrons explicitly. The fermionic effective action S^ has essentially the same 
structure as the action of the conduction electrons in Eq. (4) , considering that antiferromagnetic spin fluctuations are 
frozen to be Hi ■ f — T3 in Eq. (4), and gauge fluctuations are ignored in Eq. (15). When the spinons are condensed, 
the effective action in Eq. (15) is reduced to Eq. (4) with n, • f = t 3 . In the condensed phase gauge fluctuations 
can be ignored in the low energy limit because they are gapped due to the Andcrson-Higgs mechanism. The spinon 
action S z is equivalent to the CP 1 action of the 0(3) nonlinear a model. [5] Hence, the effective action of Eq. (15) can 
recover the insulating antifcrromagnet at half filling of the conduction electrons with Fermi-nesting. 



III. MEAN FIELD ANALYSIS 



Now we perform the saddle-point analysis to obtain possible phases of the Kondo lattice model. The mean field 
phases will be determined by condensation of the bosonic spinons. We ignore the U(l) gauge fluctuations in the mean 
field analysis. Later, these gauge excitations will be allowed beyond the mean field approximation. Then, we consider 
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the effective mean field action 

Z M F = / Dz i<7 D~lp ia . exp [- { / dT (E t E 1pl a d T 1p ilT - tE r E Iptrlpja 

icr ijcr 

icr icr icr ijcr 

+aEEi z -i 2 - 1 ) + j E^+ < E^-)}] (16) 

i <x i ij 

with A = The hopping parameters Ft, E r , F r and the effective chemical potentials A, /i can be written from 
the saddle-point equations as 

JEt = j E(i^ z -i 2 )' JF * = E^ a ^-)> 

(T CT 

^ = E(4^) ; ^ = E <^>> ■ 

CT CT 

where <5 is hole concentration. 

Integrating out fermions tl>i a and bosons z i(T in Eq. (16), we obtain the following expression for the free energy 

1 v^v^, , / iE t uj n - ii + E r et <rJ K S \ 

\ ° J « S i^ n -, + E r ei +Q ) 

+ ]3 E E ln (y n n + F r e z k + A) + E(^^t + 0£ r F r ) + E(M[1 - <S] - A). (18) 
/ccr k 

Here and ej: are the bare dispersions of chargons and spinons, respectively. to n (fi„) is the fermionic (bosonic) 
Matsubara frequency. Minimizing the free energy in Eq. (18) with respect to E t , F tl E r , F r , A, and /i, we obtain the 
self-consistent mean field equations 

to 2 JJ 



J-D t~> n _ .7 "n 



n„ j ' F r e + X' 



i, 2 ^ ibJ n [(iE t u n - /U + E r ef +Q ) + (iE t uj n - fj, + E r ef)] 



~k @ oj n ( iE tU) n - (i + E r ef)(iE t uj n — fi + E r ef +Q ) - (JkS) 2 ' 



J-D P n -f 



n n -J^n + FrC + X 



r 2-^1 R 2-^ Z7 J>\(iTP.,., .. J_ TP J> 



k 
l-D 



$ ( iE t^n - H + E r El)(iE t Lu n - /i + E r el +Q ) - (J K S) 2 ' 

1 = // eD(e) |g^ + Ve + A' 
t g = V 2 V [^ t ^ n - M + g r 6^ +Q ) + (zg t ^»-/j + £: r ef)] 

fc h ( iE t"n - V + E r et)(iE t u n - M + E r 4 +Q ) - (J K Sf ' 

Here ^ fc is replaced with J D D deD(e) in the bosonic equations, where D(e) is the density of states for the bosonic 

spectrum e|. J2 k in the fermionic equations means sum over the folded Brillouin zone. The factor 2 in the 1/(3 terms 

comes from the spin degeneracy. The chargon spectrum ef is given by the electron bare dispersion in the tight binding 

approximation, i.e., ef = —2t(cosk x + cosk y ). Furthermore, a constant density of states D(e) = 1/2D will be used, 
where D — At is a half of the band width. 
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In Eq. (19) we should introduce energy cutoff in the frequency integrals for E t and F t . Note that the usual 
momentum cut-off D is introduced in the integrals for E r and F r to prevent the divergence. Quite similarly we 
also introduce an energy cutoff J in the integrals for E t and F t because it corresponds to inverse of lattice spacing 
in the frequency space. When evaluating the frequency integrals for E t and F t , we first divide the integrals into 
two parts, divergent and divergent-free parts. We calculate the divergent parts within the energy cutoff, but for the 
divergent-free integrals we perform the Matsubara summation without the energy cutoff. 



A. Kondo insulator 



We solve the mean field equations [Eq. (19)] at half filling of the conduction electrons, where the chargon chemical 
potential \i is zero due to the particle-hole symmetry. The Fermi-nesting induces a gap corresponding to the Kondo 

coupling in the chargon spectrum. The excitation spectrum of the chargons is given by E^ — \J (E r e^) 2 + (JkS) 2 . 
Performing the Matsubara summation in the frequency integrals in Eq. (19), we obtain 



2J 1 f u , \/F r e + A , .0 F r e + X, 

JE t = h - / deD{e)^-^ — ; — coth(-J — ^ ) 

irF t J J_ D y > (f )3/2 k 2 y ^ ' 



J 

D 



2J If, VF r e + X 



nF t 2D J J_ D (S)3/2 



DE r = - ( D deD(e) i coth(^ 4 Fr€ p + A ) = (T ^ 0) (° de— 

J_ D K %^^E V 2V ^ > (T ^ 0) 2Dj_ D 4 



_ _ / F r £+\ 

J \l El V J " " J \J 11 



r V E t Et {T ^ 0) 2Dj_ D E t ^E 2 e 2 + (J K S) 2 ' 



° deD(e) p coth(^/ fre F +A ) = (M -L (° de p , (20) 



where we use a constant density of states, D(e) = 1/2D. The first terms in the equations for E t and F t come from 
the divergent parts. 

Performing the momentum integrals in Eq. (20), we obtain the following expressions for the self-consistent mean 
field equations of E t , F t , E r , F r , and A 

2J [J[X + DF r \y/ 2 - [J[X - DF r \y/ 2 

JE t — — H tjt, , 

nFt 3DJF r F? /2 

23 DE ry J D 2 E 2 + (J K S) 2 + {J K S) 2 sinh" 1 (DE r /J K S) 

J rt = — — + 



DE r = 



DF r 



irEt 2DE 2 E r 
(2A - DF r )y/J{X + DF r ) - (2A + DF r )^ J(X — DF r ) 
3DF 2 y/V t : 

DE r \J D 2 E 2 + {J K S) 2 - (J K S) 2 smh' 1 (DE r /J K S) 
2DE t E 2 . ' 



1 = y/J(X + DF r ) - VJ(X - DF r ) 
DF r JY t 



(21) 



Although it is not easy to obtain analytic expressions for the mean field parameters as a function of Jk/D and J/D, 
we can find the quantum critical point where the bosonic spinons begin to be condensed. The spinon condensation 
occurs at A c = DF TCl where c denotes "critical." Inserting this condition into Eq. (21), we find the quantum critical 
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FIG. 1: (Color online) The mean field parameters at half filling are shown as varying JkS/D with J/D — 0.01. The quantum 
critical point defined by X c /D — F rc — in (a) coincides exactly with the analytic result JkcS/D « 0.09. The effective hopping 
parameters F r and E r become zero at a certain value of the Kondo coupling beyond its critical point, which implies that our 
spin decomposition cannot cover the whole range of the Kondo lattice model. 



point defined by 

3 I 3JkcS 24J 28J J 2 J 2 4 

F tc yi + (-p-) -^^--^ + 8-, F rc -—-, E tc _-__ + _ 

1 -if3J Kc S\ 2 , u 1 \ 24 J 2 20 J 8 J 

^=3' ^( — ) (3J^)=^-V^ + 3^ (22) 

Solving Eq. (22), we find that the quantum critical point JkcS/D depends on the value of J/D. In the case of 
J/D = 0.01 we obtain JkcS/D w 0.09 > that completely coincides with the result obtained by solving Eq. (21) 
numerically as shown in Fig. 1. Decreasing J/D from this value, the critical value becomes larger. In the case of 
J/D = 0.1 we find JkcS/D —0.10 < 0, indicating that there is no phase transition at half filling, and only the 
Kondo insulating phase (see below) appears. Increasing J/D further, the critical value gets more negative. The 
condition for this quantum transition can be obtained from the boson sector of the mean field action in Eq. (16) 
without any detailed calculations. Since the boson Lagrangian coincides with the rotor model, more precisely, the 
CP 1 Lagrangian of the 0(3) nonlinear a model, one can read the transition condition (D / J)F t F r w 1 from the boson 
Lagrangian itself. Actually, it can be seen from Eq. (22) that the mean field condition for the spinon condensation 
given by A c = DF rc coincides with (D/ J)F tc F rc = 2 exactly. 

The mean field parameters in this effective Kondo action at half filling is shown in Fig. 1, where J/D = 0.01. For 
Jk < Jkc the condensation of bosonic spinons occurs, indicating the emergence of an antiferromagnetic order for the 
localized spins. On the other hand, for Jk > Jkc the spinons become gapped, implying that the localized spins are 
disordered and the antiferromagnetic order vanishes. Increasing the Kondo coupling strength, the localized spins are 
strongly affected by the conduction electrons. Thus the effective hopping parameters F t and F r decrease as increases 
Jk (Fig. 1). As F t and F r decrease further so that (D/ J)F t F r <C 1 for large Kondo couplings, quantum fluctuations 
of spinons get stronger, destructing the antiferromagnetic long range order of the localized spins. 

Meanwhile the Kondo hybridization {p,i ■ {cl^Tatr' Cio')) is nonzero in both phases because 

(4t(^ ' r)aa-Cia') = {c\ a U laa T Sa ( 3 U} pa ,C lcJ ,) = (^L^a/sV^) = (V^T ~ ^ 0, (23) 
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TABLE I: Quantum phases at half filling in the Kondo lattice model 

Jk < Jkc Jk > Jkc 

Antiferromagnetic insulator Paramagnetic insulator 

(Zia) <Z 1CT ) = 

(4 a T 3a f)Cig) + (4 a T 3a f3Cig) = 

(V>L T "3 a /3^»/3) ± {jjjsSo^M ^ 



which is the hallmark of the present strong coupling approach. In a different angle one may view this as an assumption 
in our strong coupling theory. At half filling, due to the Fermi-nesting the chargon excitations are gapped, thus both 
phases are insulators. 

When the localized spins form an antiferromagnetic order in Jk < Jkc with the condensation of bosonic spinons, 
(zur) 7^ 0, the conduction electrons also exhibit an antiferromagnetic order through the Kondo couplings with the 
localized spins. One can see this antiferromagnetic order from 

(4 a T3aa>C ilT >) = (i>l a U} aa T 3aa ,U ia , 1p i0 ) W (V>]a^M (U} ar7 T 3rTrT ,U ilT , f}) 

= (Vj T ViT - V^iX^T " z l^l) + 0, (24) 

if the easy axis anisotropy is assumed. In the easy plane limit one finds (c^Tiaa' c ia >) ^ or (c ia T2aa'Ci(r') ^ 0. 
Since the conduction electrons form an insulator due to the Fermi-nesting at half filling, the phase of the conduction 
electrons is an antiferromagnetic insulator. 

When the localized spins are in a disordered phase, the antiferromagnetic order of the conduction electrons also 
vanishes as (zi a ) = 0. As a result, for Jk > Jkc the phase of the conduction electrons is identified as a Kondo 
insulator because the conduction electrons are still gapped due to the Fermi-nesting. In the Kondo insulator the 
origin of the excitation gap is the Kondo hybridization, not the antiferromagnetic ordering. Figure 1 shows that mean 
field analysis for the effective action of the Kondo lattice model exhibits a second order phase transition from an 
antiferromagnetic insulator to a Kondo insulator, as increases the Kondo coupling strength Jk- The possible mean 
field phases at half filling are summarized in Table I. 

Although the continuous quantum transition between the two insulating phases was obtained in the mean field 
approximation, it should be considered as an artifact of the mean field analysis because instanton excitations of 
compact U(l) gauge fields[12] cause confinement of the massive spinons and chargons in the Kondo insulating phase 
beyond the mean field level. From the seminal work of Fradkin and Shenker[13] we know that there can be no 
phase transition between the Higgs and confinement phases. The order parameter discriminating the Higgs phase 
from the confinement one has not been known yet. [14] In this respect only a crossover behavior is expected. In 
this study the antiferromagnetic state corresponds to the Higgs phase because the phase is characterized by the 
spinon condensation, while the Kondo insulating state coincides with the confinement phase. Applying Fradkin and 
Shcnker's result to the present problem, we conclude that the second order phase transition turns into a crossover 
between the antiferromagnetic insulator and the Kondo insulator. This crossover picture is reasonable, considering 
that the chargon excitations are gapped in both phases. 

Note that the present spin decomposition is not allowed for all values of Jk/D because the renormalized hopping 
integrals tF r and tE r for the spinons and chargons, respectively, become zero above a certain value of Jk/D (Fig. 1) 
in the Kondo insulator. Solving the mean field equations (21) in the limit of E r — > and F r — ► 0, one can determine 
the value of JkS/D resulting in E r = and F r = from the following conditions 

These equations give J K S/D = 0.19 and E t = 8.99 for J/D = 0.01, and J K S/D = 0.13 and E t = 2.27 for J/D = 0.1. 
Both the spinon and chargon bands become flat above this Kondo coupling strength, causing these particles localized 
with (zic) — 0. We believe that this localization originates from our strong coupling approach. We interpret the 
localization as the breakdown of our spin decomposition. 



B. Heavy fermion metal 



Now we consider a hole-doped case where the Fermi-nesting is destroyed. We can expect the metallic behavior 
of chargons. Introducing the electron chemical potential, we obtain the self-consistent mean field equations for the 
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(25) 



The mean field equations in the spinon sector remains the same as those in Eq. (20). One can recover Eq. (20) for 
half filling by setting 5 = and fj, = in Eq. (25). 

In the doped case two kinds of phase transitions are expected to appear. One occurs in the spinon sector, charac- 
terized by the spinon condensation, thus associated with an antifcrro- to para- magnetic transition of the localized 
spins. The other can appear in the chargon sector, not understood by condensation of an order parameter since there 
is no order parameter in this fermion part. The phase transition is an insulator to metal transition of the chargon 
excitations, occurring when the gap in the chargon spectrum vanishes. From the last equation in Eq. (25) one can 
see how the chemical potential changes as a function of d and Jk, given by 



H = — \f (E r D) 2 5 2 + (J K S) 2 



(26) 



for S > and /j, = for S = 0. This means that as soon as holes are doped in the conduction band, the chemical 
potential that lies between the upper and lower "conduction" bands jumps to the lower band, thus metallic properties 
of the conduction electrons appear. 

The most important question in this study is how the antifcrro- to para- magnetic transition of the localized spins 
arises when the conduction electrons become metallic away from half filling. This quantum transition is driven by the 
spinon condensation. Performing the momentum integrals in the first and second equations in Eq. (25), we obtain 
the following expressions for the mean field equations 



JFt ~ ^Et + E?V DE r ) 



+ 



2DE 2 E r 



DE r y/D 2 E 2 + (J K S) 2 + /zyV 2 - (JkS) 2 
+ (J K S) 2 [sinh- 1 (DE r /J K S) - sinh-VM 2 - (J K S) 2 /J K S)] ), 

DFr = 2IjkE 2 ( DErVD2E? + {JkS? + ~ {JKS)2 

- (J K S) 2 [sinh" 1 (DE r / J K S) - smlT 1 ^ 2 - ( J K S) 2 / J K S)] ). 



(27) 



Using the equations for E t , E r , A in Eq. (21) and Eq. (27) with Eq. (26), we can find the quantum critical point 
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FIG. 2: (Color online) The mean field parameters away from half filling (6 — 0.5) are shown as varying JkS/D with J/D — 0.1. 
The quantum critical point defined by X c /D — F rc = in (a) coincides exactly with the analytic result JkcS/D w 0.36. The 
chemical potential agrees well with the analytic result fl = — \J (E r D) 2 5 2 + (JkS) 2 . Even away from half filling the effective 
hopping parameters F r and E r also vanish at a certain value of the Kondo coupling beyond its critical point. 



associated with the magnetic transition, given by 




Note that Eq. (28) is not reduced to Eq. (22) in the S — > limit because there is a chemical potential jump at half 
filling. 

Let us consider the doped case with 5 = 0.5. For J/D = 0.1 wc find Jk c S/D « 0.36, indicating the existence of 
the phase transition from an antiferromagnetic metal to a paramagnetic metal away from half filling. Remember that 
there is no phase transition at half filling when J/D = 0.1. We identify this antiferromagnetic metal to paramagnetic 
metal transition as the quantum phase transition in the Kondo lattice model. Possible mean field phases away from 
half filling are summarized in Table II, basically the same as Table I except that the phases are metallic rather than 
insulating. The mean field parameters away from half filling is shown in Fig. 2, where J/D = 0.1. Note that the 
transition point in Fig. 2, obtained by solving the mean field equations (21), (26), and (27) numerically, agrees 
completely with the analytic calculation. 

A standard way of interpreting the quantum transition uses the ground state wave function. For Jk < Ji< c , 
since (zi a ) ^ and (ipl a T3 a pipip) ^ 0, the ground state |r?i) for a localized spin at site i can be written as |fij) = 
tiaf\AF} + tiks\KS), where \AF) is the antiferromagnetic state with weight uaf, and \KS) the Kondo singlet 
state with weight tlks- We emphasize again that the Kondo hybridization still exists for Jk < Jkc- The ground 
state of the conduction electrons for Jk < Jkc is given by \ci a ) = caf\AF) + cks\KS) in the same way, where 
the antiferromagnetism and Kondo hybridization result from the Kondo interaction. Increasing the Kondo coupling 
strength, the antiferromagnetic long range order vanishes due to the Kondo hybridization, thus the ground state for 
the localized spin at site i turns into \fti) = \KS). In the same way the ground state for the conduction electrons is 
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TABLE II: Quantum phases away from half filling in the Kondo lattice model 

Jk < Jkc Jk > Jkc 

Antiferromagnetic metal Heavy fermion metal 
(z~) + (z iCT ) = 

(4 a T 3a gCip} (4 a T 3a gCig} = 

(V>L Tariff) / I^JSapjHp) 



given by \c ia ) = \KS). 

This discussion implies that fluctuations of the Kondo singlets are not critical in this transition, and only antifer- 
romagnetic spin fluctuations are critical to drive the quantum transition via the Kondo interaction. This picture is 
consistent with the single impurity problem, where the Kondo coupling always causes the Kondo singlet ground state. 
The presence of the Kondo singlets in both phases originates from the strong coupling approach, where the Kondo 
coupling term is solved first, thus allowing the Kondo singlets in both phases. This leads us to conclude that the 
volume change of Fermi surface should be continuous, which can be checked from the fact that the chemical potential 
varies continuously across the transition. 

However, there can exist another solution for the mean field equations (21), (26), and (27). If we assume E r = 
and F r = for JkS/D larger than the quantum critical point JkcS/D, these mean field equations become 

(F t E t + -)F t = l, F t Et = -, X=4r, 

V A"/ TT F t 

yielding F t = 7r/4, E t — &/n 2 , and A = 4J/n. This solution has an interesting physical interpretation although 
our spin-decomposed effective action for the Kondo lattice model is not available above the value of Jk where 
E r — > and F r — > 0. Note that this solution cannot be compatible with the spinon condensation because the spinon 
condensation occurs when {D/J)F t F r > 2 is satisfied. Thus, it can be allowed only in the paramagnetic phase. 
Approaching the quantum critical point from the antiferromagnetic phase, the effective hopping parameters remain 
finite, i.e., E rc ^= and F rc ^= while they are zero approaching the quantum critical point from the paramagnetic 
phase. There is discontinuity for E r and F r , given by AE r = E r (Jx — > Jkc — 0) — E t (Jk — * Jkc + 0) = E rc and 
AF r = F t {Jk — > Jkc — 0) — F r (Jx — > Jkc + 0) = F rc . We expect that this abrupt change in the hopping parameters 
may be related with discontinuity in the volume change of the Fermi surface[l] even if the quantum transition is the 
second order described by the spinon condensation. Although the delocalized solution (E r ^ 0, F r ^ 0) near the 
quantum critical point is expected to be a genuine solution in the mean field level, it will be meaningful to consider this 
localization solution, considering gauge fluctuations beyond the mean field approximation. Later, we will comment 
on this issue shortly. Since we don't have any clear physical picture for this localization, we do not go further based 
on this mean field solution. 



IV. BEYOND THE MEAN FIELD APPROXIMATION 



A. Effective field theory 



To examine non- Fermi liquid physics near the quantum critical point, it is necessary to obtain an effective continuum 
action of Eq. (15). It is important to introduce non-perturbative effects of Kondo interactions in the continuum action. 
We rewrite the mean field action for the fermion sector in Eq. (16), and diagonalize it as 

omf _srsr ( ,/,t ,/,t "\ ( iEtUJn - v + Er 4 aJ nS \ ( i> Un ka \ 

V 1+unk* 1-u n k* ) { iE t uj n -li- El ){ V ^ nka ) ' (2yj 



(30) 



where the rj± LUn ka fcrmions are given by the unitary transformation of the ipuj n ka fermions in the following way 

V+LO n ka \ _ ( COS?? Wnfe -(7Sini? Wnfe \ / 1p LUnk <r 
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Here = \J (E r ef) 2 + (JkS) 2 is the quasiparticlc energy obtained before, and cosi9 Wn fc, sm-d bJn k are coherence 



factors, given by cos 2 $ w „fc 



1 



and sin 2 ■& UJn k = h 



1 



Expanding the quasiparticlc band in the long wave length limit, we obtain 



E 



~~Jk~S JkS ^ x ^ ^' 



JkS 



(31) 



where the terms beyond the fourth order are ignored in the long wave length limit. Inserting the above into Eq. (29) 
and performing the Fourier transformation, one can obtain a low energy continuum action for the r\- a fcrmions 



S i>= J dr J d 2r \^2{v-a( E v\. d r - i&a T ] - \i - iA T )r\_ 

+ ^-l(v-^-aV.I 2 ) 



(32) 



where M^ 1 = (DE r ) 2 /8JkS and E t is replaced by E rj . Note that the effective mass M v of the renormalized 
conduction electrons r\- a is proportional to the Kondo coupling. The U(l) gauge field a M is introduced by shifting 
the three momentum fc M as — aa^. The empty high energy band for the rj +a fermions is ignored in the low energy 
limit. 

Performing the continuum approximation for the boson sector in Eq. (15), the resulting effective field theory is 
given by 

Seff = Sr) + S z , 

S v = J dr J d 2 r\^2(r]l(E r] [d T - iaa T ] - fi - iA T )r] a + 7^H(V - iaa - iA)j] a \ 



f\r J d 2 r[Y / (^\(d T - l a T )z a \ 2 + ^J(V- i d)z a \ 2 + rn 2 z \z a \ 2 ) 



+^(£ki 2 ) 2 +^n 2 w 2 



(33) 



where M~ x = tF r /2 and F t is replaced by F z in the spinon part. Here the "— " symbol in r)- a field is omitted for a 
simple notation. The unimodular constraint in the spinon sector is softened via their local interactions u z . The w z 
term is phcnomenologically introduced, associated with the easy-plane anisotropy. 
The above effective action can be simplified via the following scale transformation 



^/FjJ 1 

Performing the field-transformation for the boson sector accordingly. 

a 



^2M z r. 



a' T = ^F z /Ja T , 



a/4 



a = 



^/2M~ Z 



V2M~z 



the effective spinon action is given by 



S z = £ dr' J d 2 r'[^(\(d^- i a' T )z' a \ 2 + \(V , -za')z'^ + m z 2 K 



'2| J \2 



(34) 



(35) 



+^(El^l a ) a + ^-KlVil 



(36) 



with 



VK/j 



m' z 2 =m 2 z , u' z 



2M Z 



2M Z 



(F z /jy/ 2 



UK 



(F z /jy/ 2 
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The effective chargon action can also be obtained as 



9 dr' f d 2 r' [j2(ri'J(E' v [d T , - iaa' T ] - // - iA' T W a + ^|(V' - iaa' - il'WJ 2 



with the scale transformation for the fermion sector 



_ E v 



M' 



Mr, 

2M Z 



A' 



(F z /J) 



1/4 



^/2M Z 



n' = fj,, A' T = A T . 
As a result, we find the effective field theory 



S ef f = J dr J d 2 r\^2(ril(E v [d T - iaa T \ - ^ - iA T )r] a + - iaa - iA)r) a \ 2 ^ 
+ fdrj d\ [J2(\(dr - l a T )z a \ 2 + |(V - ia)z a \ 2 + m\\z a f) + y (£ ^ + 



wM\ \ z i\ 



(37) 



(38) 



(39) 



where the prime symbol is omitted for a simple notation. In the gauge action the former with an internal gauge 
charge g of the r\ a and z a particles is the Maxwell term resulting from high energy fluctuations of the r\ a and z a 
particles. [15, 16] The latter with the quasiparticle " renormalization" E n and the Fermi velocity of the chargons 
is the Landau damping term representing dissipativc dynamics of gauge fluctuations, which come from particle-hole 
excitations of the rj a fermions near the Fermi surface. [15] Since the time component of the gauge field mediates local 
interactions due to the T) a polarization[15] which are irrelevant in the renormalization group sense, it can be ignored 
in low energy limit. 



B. Antiferromagnetic metal 

For weak Kondo couplings Jk < Jr c (ml < 0) the bosonic spinons become condensed, leading to an antiferromag- 
netic order of the localized spins. An antiferromagnetic metal appears for the conduction electrons away from half 
filling. Gauge fluctuations are massive due to the Anderson-Higgs mechanism, thus safely ignored in the low energy 
limit. Although the gauge fluctuations are irrelevant in the renormalization group sense, they play an important role 
in confining the fermionic chargons with the bosonic spinons to make usual conduction electrons. This can be seen 
from the unitary gauge. If the easy plane limit w z < is considered in Eq. (39), the spinons can be treated as 
z a = (l/v / 2)e l< ^' T . The unitary gauge means a M = a M — causing an excitation gap for the a M fields. In this 

unitary gauge the phase degrees of freedom appear in the chargon sector, and these phase fields can be gauged away 
from the gauge transformation c a = e~ %a ^^ r\ a . Thus, low energy excitations are antiferromagnons e J ^ T_< ^^ in the 
localized spins and electron excitations c a in the renormalized conduction band instead of the fermionic chargons r\ a . 
These conduction electrons feel week staggered magnetic fields due to the antiferromagnetic ordering of the localized 
spins. 

C. Quantum critical point 

As increases the Kondo coupling, the antiferromagnetic metal approaches the quantum critical point where the 
antiferromagnetic order vanishes. Critical boson fluctuations renormalize gauge dynamics[17] in the critical field 
theory 

S c = jT dr J d 2 r\j2(vt(E n d T - n - iA T )r) a + ^H(V - iod - iA) V<J \ 2 ^ 
. 1 (E n \w n \ N z w qiqj\ 
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where N z = 2 is the flavor number of the spinons. Integration for the critical spinons should be understood in the 
renormalization group sense. Since critical boson fluctuations yield a term linearly proportional to momentum q in 
the gauge action, the dynamical critical exponent is obtained as z — 2. [18] 
In the random phase approximation[19, 20] the free energy is given by 
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(41) 



where D(q,uj) = {^—i^f 1 ^ + 1S the gauge kernel in the real frequency uj, and uj c is an energy cutoff. The specific 

heat is obtained to be 



1 E„ 



\dT 2 )v n 2 v„N z v v N z ir 2 v v N z 8tt 2 T T 



1 T T 
TlnT = -— - — In — 



(42) 



with an energy scale T = ) • Thus, the specific heat coefficient 7 has a singular dependence 7 = Cy/T cx 

— InT in the T — > limit. The logarithmic divergence also can be seen in the two dimensional itinerant antiferro- 
magnet, where its critical field theory is characterized by the dynamical exponent z = 2. [20] 
In the one-loop level the imaginary part of the fermion self-energy is given by 



V(MZ)= I dco J ^p[n(w) + !][!-«)] 
) a/J -^M)lmI>( g ,a;)<5( e 2- 



(k + k') a (k + k')t 
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M 2 N z €k ' 



(43) 



where = k 2 /2M v is the energy dispersion of the r] a fermions, and N v is the density of states at the Fermi energy. The 
scattering rate can be obtained from the self-energy expression with an additional q 2 factor in the integrand. [15] Since 
the imaginary part of the fermion self-energy is linearly proportional to the fermion dispersion, the dc conductivity[21] 
is given by a c <~ T~ 2 . 



D. Region of strong Kondo couplings 

Increasing Jk further from the quantum critical point, an anomalous metallic phase appears. Integrating out the 
gapped z a excitations in Eq. (39), we obtain the Maxwell term for gauge fluctuations 

Snfl = j£ dr J d 2 r^2(r]UE v d T - ^ - iA^ + - icra - iA)r] a \ 2 ^ 

l\r^/9 2 . E, q \uj n \\( qiqj\ 

where the dynamical critical exponent is z = 3. The effective field theory of Eq. (44) is well known to cause non- 
Fermi liquid physics due to scattering with massless gauge fluctuations. The imaginary part of the fermion self-energy 
is given by uj 2 ^ at the Fermi surface, implying that its real part also has the same frequency dependence via the 
Kramer's Kronig relation, thus giving rise to a non- Fermi liquid behavior. [22] The coefficient 7 of the specific heat is 
proportional to — InT in three spatial dimensions and T~ 1//3 in two dimensions. [8] The dc conductivity is proportional 
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TABLE III: Physical response in the Fermi liquid to non-Fermi liquid transition 



Jk < Jkc Jk ~ Jkc Jk > Jkc 

Antiferromagnetic Fermi liquid Quantum critical point Paramagnetic non-Fermi liquid 



7 = C v /T 




to T~ 5 / 3 in three dimensions and T~ 4 / 3 in two dimensions. [15] Note that the dynamical exponent z changes from 
z = 2 at the quantum critical point to z — 3 in the non-Fermi liquid phase. Physical responses in the Fermi liquid to 
non-Fermi liquid transition for two spatial dimensions are summarized in Table III. 



When the spinon excitations are gapped, they can be ignored in the low energy limit. Thus, if the gauge fluctuations 
are suppressed in Eq. (44), Fermi liquid physics can be obtained. As the Kondo coupling constant increases, the 
effective chargon mass M v = 8JkS/ (DE r ) 2 becomes heavier and gauge fluctuations are suppressed because l/v n ~ M v 
in Eq. (44). This may give rise to the Fermi liquid physics in the case of large Kondo couplings. In this scenario 
the non-Fermi liquid is expected to turn into the Fermi liquid continuously. In our mean field analysis M z — > oo and 
M n — ► oo were found at the point where E r — and F r = 0. Hence, the spin decomposition scheme cannot cover the 
whole range of the phase diagram so that we were not able to recover the Fermi liquid phase. 

There is another possibility associated with the confinement-deconfinement transition due to the compactness of 
the U(l) gauge field in the present problem. Note that we did not take into account instanton excitations in the 
previous discussion. In two space and one time dimensions there is no deconfined phase owing to proliferation of 
instanton excitations when only gapped fermion or boson excitations exist. [12] However, the presence of gapless 
matter fields was recently argued to allow a deconfined phase. [16, 23-28] Non- Fermi liquid phase corresponds to 
the deconfined phase which gapless fermion (rj a ) excitations make stable against instanton excitations. [16, 23] The 
present quantum critical point is identified as the deconfined quantum critical point [24] that can be stable due to 
critical boson (z a ) excitations [25] and gapless fermion excitations. [16, 23] On the other hand, Fermi liquid corresponds 
to the confinement phase. Instanton condensation leads to confinement between the r\ a fermion and the z a boson 
to make an electron c CT = V aa ir\ a i. We don't know whether the non- Fermi liquid phase is stable against instanton 
excitations or not. If it is stable, the confinement-dcconfincmcnt transition corresponding to the non-Fermi liquid 
to Fermi liquid transition would occur in the strong Kondo coupling region. The nature of this transition may be 
KT-likc (Kostcrlitz-Thouless).[29] On the contrary, if the non- Fermi liquid phase is unstable against the confinement, 
the parameter region of the non-Fermi liquid phase would shrink to vanish. Then, the quantum critical point will 
coincide with the point where localization occurs with M z — > oo and M v — ► oo. 



Our approach has some analogies with that of Ref. 30 where bosonic spinons arc used for the localized spins, 
resulting in charged fermions for the Kondo resonances. However, there are several important differences between our 
approach and that of Ref. 30. Ref. 30 starts from the Kondo-Heisenberg lattice model 



E. How to recover the Fermi liquid phase 



V. 



DISCUSSION AND PERSPECTIVES 




(ij)aa' 



(45) 



~* it-* 

where Cka represents a conduction electron with momentum k and spin a, and Si = ^&] cr T (7(7 '&j cr > is the boson 
representation of the localized spin Si. Performing the HS transformation for the particle- hole channel in the Kondo 
coupling term and the particle-particle channel in the Heisenberg interaction term, Eq. (45) reads 





(46) 
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FIG. 3: (Color online) A schematic phase diagram of the Kondo lattice model based on our effective theory 

where the on-site bond variable x% 1S a Grassman field associated with the Kondo resonance and the bond variable 
Aij is introduced to keep short range antiferromagnetic correlations. 

The crucial difference between two approaches lies in the HS decoupling scheme of the Kondo interaction term; 
Rcf. 30 allows three kinds of matter fields that correspond to two fermions c^a, Xi an d one boson 6 icr while our 
decomposition introduces only two kinds of matter fields, one fermion %pi a and one boson Zi a . In fact, Xi fermion 
corresponds to ip^ fermion while ipil fermion is not allowed in Ref. 30. Since Xi field follows fermion statistics and 
the condensation of fermions is not possible, the conventional mean field analysis in the slave-boson approach[3] is not 
applicable. Recently, there has been a progress in this spin-boson approach for the single impurity problem based on 
the non-crossing approximation scheme of the U(l) slave-boson theory,[3] although its extension to the Kondo lattice 
model has not been reported yet. [31] 

A schematic phase diagram based on the effective Kondo action in Eq. (15) is shown in Fig. 3, where horizontal 
axis is the Kondo coupling strength and vertical is temperature. "AF" represents the antiferromagnetic metal and 
"HF" the heavy fermion metal. The dashed line denoted by Tk represents the Kondo temperature for the Kondo 
singlets to form, where (atp^ipia) ^ below Tk while {aipj^tpia) = above Tk- The solid line shows the second order 
antiferromagnetic transition related with the spinon condensation in the present theory. This phase diagram is quite 
similar to that of the HMM theory since the Kondo hybridization always exists in both the antiferromagnetic and 
heavy fermion phases below the Kondo temperature, and antiferromagnetic ordering is associated with the quantum 
transition. However, an important difference between these two theories can be found near the antiferromagnetic 
quantum critical point. 

In the HMM theory low energy elementary excitations at the quantum critical point are critical antiferromagnetic 
fluctuations with spin quantum number 1, described by 

Shmm = ^ y^(r|cj w | + q 2 )n(q,w n ) ■ n(-q, -UJ n ) + V(\n\), (47) 

where the damping term with a damping coefficient T comes from gapless electron excitations near the Fermi surface 
and V(|n|) = Jq dr J d 2 r 1 rj L \ri\ 4: + •■■ is an effective potential for the spin-fluctuation order parameter n.[2] On the 
other hand, at the quantum critical point of the present approach the critical antiferromagnetic fluctuations are 
fractionalized into critical spinon excitations with spin quantum number 1/2 due to strong Kondo interactions. The 
critical field theory is given by 

Sdqcp= £ dr J d 2 r[Y,(\(dr - ia T )z a \ 2 + |(V - m)z CT | 2 ) 

+y(E w 2 ) 2 + ^M 2 hi 2 ] + \ E + f q ) i 5 * - q ~fh a » (48) 

where the gapless renormalized fermions r\ a are integrated out, causing dissipation in gauge dynamics. Here "DQCP" 
means the deconfined quantum critical point [24, 25] discussed in the previous section. 
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It is interesting to see that both critical theories in Eqs. (47) and (48) are characterized by the same dynamic 
critical exponent z = 2. Thus, two critical theories show similar critical physics although the low lying excitations 
are completely different. Actually, the HMM theory can also explain the logarithmic divergence in the specific heat 
coefficient in two spatial dimensions. [1] Since two spatial dimensions lie in the upper critical dimension due to the 
dynamic critical exponent, universal scaling for the spin susceptibility does not appear in both critical theories. Thus, 
it is impossible to compare an anomalous critical exponent in one theory with that of the other. However, there exists 
one crucial difference; if we have the deconfined quantum critical point, the non- Fermi liquid metal as a deconfincd 
critical phase can appear above the deconfined quantum critical point ( Jk > Jkc) while this non-Fermi liquid phase 
cannot be allowed in the HMM theory. According to our effective field theory, it is clear that there should be a 
crossover between the quantum critical region and non- Fermi liquid phase as decreases temperature in Jx > Jkc- 
The crossover temperature T* depends on the spin gap \ — DF r as T* ~ X — DF r . This crossover should appear in the 
upturn behavior of the specific heat coefficient from 7 ~ — In T to 7 <~ T -1 / 3 since the effective field theory changes 
from z — 2 to z = 3 during the crossover. The upturn behavior was also discussed in Ref. 30, but the mechanism is 
different. 

Let us discuss dissipation in both the weak coupling critical theory of Eq. (47) and the strong coupling one 
of Eq. (48). The strong coupling theory may be derived from the weak coupling one using the CP 1 representation. 
Remember that the 0(3) nonlinear a model can be mapped onto the U(l) gauge theory within the CP 1 representation, 
as discussed before. The main point of this CP 1 decomposition is how dissipative dynamics of spin fluctuations n(r, r) 
in the HMM theory [Eq. (47)] is transferred into that of gauge fluctuations a(r, r) in the CP 1 gauge theory [Eq. (48)]. 
In the context of the standard weak coupling theory order parameter fluctuations directly couple to gapless fermion 
excitations near the Fermi surface. As a result, dissipation effects in order parameter fluctuations appear in the 
kinetic energy term. In our strong coupling approach fractionalized order parameter fluctuations do not couple to the 
gapless fermion excitations directly. Instead, their couplings are realized indirectly via gauge fluctuations. Dissipative 
dynamics of fractionalized order parameter fluctuations is induced by damped gauge fluctuations that result from the 
gapless fermion excitations near the Fermi surface. The damped gauge fluctuations play an important role in quantum 
critical physics. If fermion excitations are gapped and the damping effects in gauge fluctuations are not taken into 
account in Eq. (48), the quantum phase transition belongs to the inverted XY (IXY) universality class in the case of 
large flavors of fractionalized boson excitations. [17] However, the presence of dissipation in gauge excitations due to 
gapless fermions changes the IXY universality class completely. Since the dissipation results in the dynamic critical 
exponent z = 2, the spacial dimension d = 2 becomes the upper critical dimension and the nature of the quantum 
transition would be a mean field-like type with logarithmic corrections. 

Inserting the CP 1 representation Ti — -^z^Tfjfj' z a i into Eq. (47), [32] Eq. (47) can be written in terms of the 
bosonic spinons interacting with gauge fluctuations. Unfortunately, the \u n \ linear (damping) term prevents obtaining 
a complete expression. Performing the HS transformation for the damping term, dissipation in order parameter 
fluctuations would appear in gauge fluctuations although the dissipative gauge action is not given by that in Eq. (48). 
This implies that damping effects due to gapless fermions are imposed in a different way for the weak and strong 
coupling theories. 

It is valuable to apply our spin decomposition to the one dimensional Kondo lattice model in order to confirm that 
the non-Fermi liquid metal with spin gap can be allowed. The effective field theory in Eq. (39) will be applicable to 
the one dimensional case. An important difference from the two dimensional case is that dissipative dynamics in gauge 
fluctuations does not appear because the gapless conduction fermions are described by massless Dirac fermions near the 
Fermi points. In one dimensional effective field theory strong quantum fluctuations coming from low dimensionality 
do not allow the spinon condensation. Furthermore, the spinon excitations are gapped because the Berry phase 
contribution disappears due to the Kondo coupling. The most crucial point in the one dimensional effective theory 
is that the gapless Dirac fermions make gauge fluctuations massive, [33] which can be seen using the bosonization 
technique. As a result, the gapped spinon excitations are deconfined [34] although the mechanism is different from the 
two the dimensional case. Moreover, the gapless fermion excitations exhibit strong superconducting correlations as the 
two dimensional case. Our spin-gauge theory in Eq. (39) allows superconducting instability because the spin-gauge 
fields mediate attractive interactions between rft and T]i fermions. The <r symbol in the gauge coupling shows that 
the gauge charge of the jjf fermion is opposite to that of the rji fermion. In this respect the non-Fermi liquid phase 
with spin gap is the two dimensional analogue of the one dimensional spin-gapped phase. 

If SU(2) gauge fluctuations are taken into account in Eq. (11) instead of U(l) gauge fluctuations, the non-Abelian 
nature of SU(2) gauge fluctuations may not allow the deconfined quantum criticality and non- Fermi liquid phase. [12] 
However, there is no consensus for the confinement problem in the context of the SU(2) gauge theory, as far as we 
know. If the deconfined non-Fermi liquid phase turns into a confinement state, spinons should be confined with 
chargons via SU(2) gauge fluctuations. Instead, electron excitations are allowed and the resulting phase may be the 
Fermi liquid. The spinons and chargons are not meaningful objects in the low energy limit. However, the spinon 
excitations may emerge as broad spin spectrum (particle-hole continuum) at high energies beyond multi-paramagnon 
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scattering according to the asymptotic freedom of the SU(2) gauge theory. [12]. 

In this SU(2) gauge theoretic description the quantum critical point would lie between the Higgs phase (antiferro- 
magnetism) and the confinement one (Fermi liquid), while it is between the Higgs phase and the deconfinement one 
(non-Fermi liquid) in the present U(l) gauge theory. It was argued that there is no phase transition between the Higgs 
and confinement phases and they are smoothly connected. [13] Then, the spin-decomposition method in the context 
of the SU(2) gauge theory cannot describe the quantum phase transition of the Kondo lattice model. In this case 
another order parameter should be considered to study the quantum phase transition of the Kondo lattice model, for 
example, the hybridization order parameter in the context of the slave-boson theory. [8, 9] 

In summary, we investigated the quantum phase transition from an antiferromagnetic metal to a heavy fermion metal 
in the Kondo lattice model. First, we diagonalizcd the Kondo coupling term in the strong coupling approach. Then, 
we derived the effective Kondo action [Eq. (15)] and performed the mean field analysis [Eq. (19)] to obtain the mean 
field phase diagram, showing the quantum phase transition from the antiferromagnetic metal to the heavy fermion 
metal. The Kondo term is always relevant so that the Kondo hybridization persists even in the antiferromagnetic 
metal, which means that fluctuations of the Kondo singlets are not critical in the phase transition. The volume change 
of Fermi surface thus is expected to be continuous across the transition. We found that softening of antiferromagnetic 
spin fluctuations leads to the quantum transition driven by the Kondo interaction in the strong coupling approach. 
Beyond the mean field level we derived the effective U(l) gauge theory [Eq. (39)] in terms of the rcnormalizcd 
conduction electrons r\ a and the spin-fractionalized excitations z a interacting via the U(l) spin-gauge fields a M . Our 
critical field theory characterized by the critical exponent z = 2 can explain the non-Fermi liquid physics such as 
7 ~ — InT near the quantum critical point. Furthermore, we showed that if our scenario is applicable, there can 
exist a narrow region of the non- Fermi liquid phase with spin gap near the quantum critical point. We also discussed 
how the present theory can recover the Fermi liquid phase, but this issue should be clarified near future. Lastly, we 
commented on the superconducting instability near the quantum critical point. This interesting possibility remains 
as a future study. 
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